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ABSTRACT 

We have observed a dramatic change in the spectrum of the formerly heavily ab- 
sorbed 'overlapping-trough' iron low-ionization broad absorption line (FeLoBAL) 
Jjj ' quasar FBQS J1408+3054. Over a time span of between 0.6 to 5 rest-frame years, 

c/2 , the Mgll trough outflowing at 12,000 km s^ 1 decreased in equivalent width by a fac- 

tor of two and the Fen troughs at the same velocity disappeared. The most likely 
explanation for the variability is that a structure in the BAL outflow moved out of 
. our line of sight to the ultraviolet continuum emitting region of the quasar's accretion 

^ ' disk. Given the size of that region, this structure must have a transverse velocity of 

between 2600 km s _1 and 22,000 km s _1 . In the context of a simple outflow model, 
CN| i we show that this BAL structure is located between approximately 5800 and 46,000 

Schwarzschild radii from the black hole. That distance corresponds to 1.7 to 14 pc, 11 
to 88 times farther from the black hole than the H/3 broad-line region. The high veloci- 
ties and the parsec-scale distance for at least this one FeLoBAL outflow mean that not 
all FeLoBAL outflows can be associated with galaxy-scale outflows in ultraluminous in- 
frared galaxies transitioning to unobscured quasars. The change of FBQS J1408+3054 
from an FeLoBAL to a LoBAL quasar also means that if (some) FeLoBAL quasars 
have multiwavelength properties which distinguish them from HiBAL quasars, then 
some LoBAL quasars will share those properties. Finally, we extend previous work on 
' how multiple-epoch spectroscopy of BAL and non-BAL quasars can be used to con- 

5-j \ strain the average lifetime of BAL episodes (currently >60 rest-frame years at 90% 

confidence). 

Key words: galaxies: active - quasars: general - quasars: absorption lines - quasars: 
individual: FBQS J1408+3054. 



1 INTRODUCTION 

Broad Absorption Line (BAL) quasars are a subset of 
quasars which exhibit strong absorption troughs from gas 
flowing outwards with velocities of thousands to tens of 
thousands of km s _1 (e.g., Allen et al. 2010). The traditional 
definition is that BAL quasars are those which have C IV ab- 
sorption which dips at least 10% below the continuum level 
for a contiguous stretch of at least 2000 km s _1 at blueshifts 
of ^3000 km s _1 relative to the quasar (Weymann et al. 
1991). Other proposed definitions (Hall et al. 2002; Trump 
et al. 2006; Gibson et al. 2009) relax some of the condi- 
tions of the traditional definition to include more cases of 
potential intrinsic absorption with the tradeoff of increased 



contamination by blended narrow absorbers (Knigge et al. 
2008) whose origin(s), especially at low outflow velocities, 
may differ from those of BAL outflows. 

BAL quasars are often classified into three subtypes: 
high-ionization (HiBAL), low-ionization (LoBAL) and iron 
LoBAL (FeLoBAL). HiBAL quasars show absorption from 
only relatively high-ionization species such as C IV and N v. 
LoBAL quasars show high-ionization features plus absorp- 
tion from low-ionization species such as Mg n and Al ill. 'Iron 
LoBAL' (Becker et al. 1997) or 'FeLoBAL' (Becker et al. 
2000) quasars show high- and low-ionization troughs plus 
absorption from excited states of Fell, Fein, or both. There 
is no commonly accepted definition for how broad or deep 
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the Mgn or Fen absorption has to be for a quasar to be 
classified as a LoBAL or FeLoBAL quasar (Trump et al. 
2006; Gibson et al. 2009; Zhang et al. 2010). However, such 
troughs should at least be broader than intervening Mgn 
absorption systems, the vast majority of which have widths 
<500 km s" 1 (Mshar et al. 2007). A minimum depth of 10% 
is appropriate for LoBAL and FeLoBAL troughs just as for 
HiBAL troughs, though in all cases this is not a physical 
limit but reflects the limitations imposed by spectral signal- 
to-noise ratios and our inability to model quasar continua 
exactly. 

BAL troughs are often highly saturated, so that the 
trough depth is determined less by the gas column density 
and more by the covering fraction; that is, the fraction of the 
continuum source occulted by the absorbing structure (e.g., 
Arav et al. 1999). Variability in BAL trough strengths is 
relatively common (see the exhaustive reference list in Hall 
et al. 2007, and Gibson et al. 2008) but large (>50%) changes 
in absorption equivalent width are rare (Crenshaw et al. 
2000; Lundgren et al. 2007; Hamann et al. 2009; Leighly et 
al. 2009; Gibson et al. 2010). 

The 2=0.848 quasar FBQS J140806.2+305448 (here- 
after FBQS J1408+3054) was first reported by Gregg et 
al. (1996), where it was misclassified as a star, and was 
first identified as a LoBAL quasar in White et al. (2000). 
Hall et al. (2002) described it as a member of the class 
of 'overlapping-trough' FeLoBAL quasars because its Fell 
absorption troughs are broad enough to overlap and yield 
almost no continuum windows below 2800 A. However, the 
covering fraction of its Fell absorption was only ~75% at 
its discovery epoch, unlike some overlapping-trough objects 
which have > 95% covering fractions. FBQS J1408+3054 
also likely has an FR II radio morphology (Gregg et al. 
2006). 

In this work we present spectroscopy of FBQS 
J1408+3054 spanning nearly 8 years in its rest frame. 
The spectroscopic dataset shows mild variability followed 
by a dramatic decline in Fell and Mgn absorption 
strength (§|5J. Substantial trough variability has previously 
been observed in the overlapping-trough FeLoBAL quasar 
SDSS J043742.81-004517.6 (Hall et al. 2002), but FBQS 
J1408+3054 is the first known case of an FeLoBAL spec- 
trum changing to a LoBAL spectrumQ We estimate a black 
hole mass for FBQS J1408+3054 in §0 In §0]we constrain 
the BAL structure's size and transverse velocity and present 
a simple outflow model which we use to constrain the BAL 
structure's location and age. We also discuss the implica- 
tions of our observations for the average episodic lifetime of 
BAL troughs in quasar spectra and for the link(s) between 
BAL subtypes and between BAL troughs and FR II radio 
morphologies. We end with our conclusions in § [5] 



We note that the absorption system D+E in NGC 4151, while 
too narrow (FWHM = 435 km s — 1 ) and low-velocity (491 km s" 1 ) 
to qualify as a BAL trough, exhibited only Mgll absorption 
in 1994-1996 but also exhibited strong Fell absorption in 1999 
(Kraemer et al. 2001) and weak Fell absorption in 2002 (Krae- 
mer et al. 2006), indicative of transverse motion of an intrinsic 
quasar absorption system such as that reported herein. 



2 DATA 

2.1 Optical Spectroscopy 

Spectroscopy was obtained using numerous telescope and 
instrument combinations (Table 1), including the Echellette 
Spectrograph and Imager (ESI; Sheinis et al. 2002) at the 
Keck telescope, the Kast Spectrograph on the Lick Obser- 
vatory 3m Shane telescope (Miller & Stone 1993), and the 
Low Resolution Spectrograph (LRS; Hill et al. 1998) on 
the Hobby-Eberly Telescope (HET). Data were obtained, 
reduced and calibrated using standard procedures (Gregg 
et al. 1996; White et al. 2000; Gregg et al. 2002), except as 
we now discuss with regard to the HET spectrum. 

We obtained a spectrum of FBQS J 1408+3054 with the 
LRS on the HET with the g2 grism (useful observed wave- 
length coverage 4335-7155 A) and a 1" slit. We obtained 
two spectrograms with a total exposure of 1510 seconds in 
2" seeing at an average airmass of 1.26. A flux calibration 
star was not observed. We used the known throughput of 
the instrumental configuration to remove the instrumental 
response before correcting for interstellar and atmospheric 
extinction. When compared to the other normalized spectra 
shown in Figure 1, the normalized HET spectrum exhibited 
a large-scale curvature resulting in a downturn in flux over 
several hundred observed A at both ends of the spectrum. 
We believe this curvature resulted from slit losses due to dif- 
ferential refraction. For the slit to include the quasar and a 
star bright enough to guide on, the long axis of the slit had to 
be oriented about 82° from the parallactic angle, resulting in 
differential refraction of ~0.75" over our wavelength range 
(Table I of Filippenko 1982, hereafter F82). Inspection of 
the results of the two individual exposures obtained showed 
that one exposure yielded a more prominent downturn in 
its spectrum and also showed more flux from a third object 
slightly offset from the centre of the slit, suggesting that 
the quasar was not well centred in the slit in that exposure. 
We therefore dropped that exposure from consideration. To 
correct for slit losses in the remaining exposure, we assumed 
that the target was centred in the slit at a wavelength of 
6500 A. We computed the offset of the target along the nar- 
row axis of the slit from 4000-7500 A using Table I of F82. 
We then used equation 8 of F82 to calculate the fraction 
of the flux captured by the slit as a function of wavelength 
and corrected the spectrum for those wavelength-dependent 
losses. 

The panels of Figure Q] show all of our spectra of FBQS 
J 1408+3054 taken between 1995 and 2009, normalized to an 
average of unity between observed 6000 A and 7000 A. The 
dashed vertical lines show the wavelengths of H/3 and Mg n 
emission in the quasar rest frame. The dotted vertical lines 
show absorption from Fen A2632 and Mgn AA2796,2803 at 
the redshift of maximum absorption (z=0.775, a blueshift 
of 12,000 km s _1 fr om the quasar rest frame). Blueward of 
Fen A2632, the spectra suffer overlapping absorption from 
hundreds of transitions of Fen. A strong decrease in absorp- 
tion strength occurred over the 955 rest-frame days between 
the two Keck epochs (MJDs 51612, in red, and 53377, in 
cyan). The absorption was still dramatically weaker in the 
Sloan Digital Sky Survey (SDSS; York et al. 2000) spectrum 
(in blue) taken on MJD 53795, 123 rest-frame days after the 
second Keck epoch, and made available in the SDSS Seventh 
Data Release (DR7; Abazajian et al. 2009). The corrected 
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Figure 1. Spectra and error arrays of FBQS J1408+3054 from 1995 (top left) to 2009 (bottom right), each smoothed with a seven-pixel 
boxcar, colour-coded, and labelled with the MJD of observation. For comparison, in all panels the MJD 51612 Keck ESI spectrum is 
plotted in red. The vertical axes are normalized F\\ all spectra have been normalized to an average of unity in the observed wavelength 
range 6000-7000 A. The lower horizontal axes are observed wavelengths and the upper are rest-frame wavelengths at z em =0.848, both in 
A. Dashed vertical lines show the wavelengths of emission from H/3 and Mg II. Dotted vertical lines show the wavelengths of absorption 
from Fell A2632 (the longest-wavelength component of the UV1 multiplet) and Mg II AA2796,2803 at z a b s =0.775, the central redshift 
of the outflow. In the left-hand panels, only minor variations in the absorption are seen to occur. In the right-hand panels, a strong 
weakening in the absorption is evident between the Keck ESI spectra taken on MJD 51612 (red) and on 53377 (cyan). This weakening 
continues in the MJD 53795 SDSS spectrum (blue) and in the MJD 55028 HET LRS spectrum (green). 
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Figure 2. Rest-frame near-ultraviolet spectra and error arrays of FBQS J1408+3054 during the disappearance of its Fell troughs. The 
spectra have been smoothed with a seven-pixel boxcar, colour-coded as in Figure [l] and labelled with the MJD of observation. The lower 
horizontal axis shows wavelengths in the rest frame of peak absorption (z=0.775) and the upper shows rest-frame wavelengths at the 
quasar z em =0.848, both in A. The vertical dashed and dotted lines are the same as in Figure [T] The dark grey shaded region shows 
the wavelength range of low- velocity Mg II absorption (squares in Figure [HJl . The light grey shaded regions show the wavelength ranges 
absorbed by Fell multiplets UV1 and UV2+UV3 at z=0.775 and by a blend of Fen UV62 and high- velocity Mgn absorption. Regions of 
absorption by Fell multiplets Opt6+Opt7 and Opt8 are also labelled, as are broad emission lines of Mgn and Fell. Note that the MJD 
55028 spectrum shows some weakening of the Fe II Opt8 emission at ~2950 A in the quasar rest frame. 
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HET spectrum, taken 668 rest-frame days after the SDSS 
epoch and shown in green in the lower right panel of Fig- 
ure!]] reveals Mg II absorption even weaker than in the SDSS 
epoch and a complete lack of Fe II UV1 absorption. 

Figure [2] displays our near-ultraviolet spectra of FBQS 
J 1408+3054 from the period of weakening absorption. Wave- 
lengths (in A) are given in the rest frame of peak absorption 
(jz=0.775) along the bottom and in the quasar rest frame 
(zem=0.848) along the top. As in Figure [T] the vertical 
dashed line shows the wavelength of Mg n in the quasar rest 
frame defined by the H/3 line (in the more recent spectra, 
the Mgn peak is blueshifted relative to the H/3 rest frame) 
and the rightmost vertical dotted line shows the wavelength 
of Mgn at the velocity of peak absorption, while the left- 
most dotted line shows absorption at the same velocity from 
Fell A2632, the longest-wavelength component of the Fell 
UV1 multiplet (Moore 1950). The left two light grey shaded 
regions show the wavelength ranges absorbed by Fen mul- 
tiplets UV2+UV3 and UV1 at 2=0.775; as expected, Fen 
A2632 traces the long-wavelength edge of the UV1 absorp- 
tion trough. Absorption outside the shaded regions arises 
predominantly from hundreds of excited Fell transitions 
whose troughs overlap in velocity space. 

The dark grey shaded region shows the wave- 
length range spanned by Mgn absorption at velocities 
i;=8700 km s" 1 to v=14,200 km s" 1 . The rightmost light 
grey shaded region shows the wavelength ranges absorbed 
by a blend of Fen UV62 + higher- velocity Mgn absorp- 
tion. This high- velocity wing of the Mgn trough extends to 
2623 A in the quasar rest frame (v=19,500 km s ), as seen 
in the HET spectrum from MJD 55028 (green). The ab- 
sorption at those wavelengths in the HET spectrum cannot 
be Fen UV62 because the Fen UV62 multiplet is intrinsi- 
cally weaker than the UV1 multiplet and the latter is not 
present in the HET spectrum. Without detailed modeling, 
it is impossible to be sure that the Fe n absorption spanned 
the same velocity range as Mgn, but the shape of the MJD 
51612 spectrum is consistent with Fe n UV62 absorption also 
extending to 19,500 km s" 1 . 

Figure |3] shows the Mgn and blended Mg n+Fe n UV62 
rest-frame equivalent widths W r for all spectra in Figure [T] 
The equivalent widths were measured relative to a contin- 
uum established using the normalized SDSS spectrum shown 
in Figure 1. In that spectrum, a linear continuum in F\ vs. 
A was drawn between the average fluxes at rest-frame 2550- 
2600 A and 3765-3820 A. In all spectra, Mgn W r values 
were measured between 2670-2720 A and Mg n+Fe n UV62 
W r values were measured between 2585-2720 A. While there 
is some uncertainty in the actual continuum level and thus 
in the measured W r , the changes in W r in Figure[3]are much 
larger than those uncertainties. As the Mg n+Fell UV62 
blend consists only of high- velocity Mg n in the HET epoch, 
we conservatively estimate that the W r of Mgn had de- 
creased by a factor of two by the HET epoch. 

2.2 Multiwavelength Photometry 

In Table 2, we present time-ordered observed and syn- 
thesized photometric magnitudes and colours of FBQS 
J 1408+3054. The photographic magnitudes listed are the 
average of the USNO-B1.0 (Monet et al. 2003) and GSC2.3.2 
(Lasker et al. 2008) magnitudes whenever available and are 
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Figure 3. Rest-frame equivalent widths in A of Mg n+Fe II UV62 
(crosses) and low- velocity Mgll (squares). Error bars shown for 
each point do not include the uncertainty in continuum place- 
ment. Around rest-frame day 400, the absorption in both species 
increased in strength quickly (over 19.5 rest-frame days) and then 
declined gradually (over 13.5 days). Between rest-frame days 948 
and 1903, the Mg Il+Fe II UV62 W r fell by nearly half. (The points 
without error bars at day 1759 are photometric estimates; see 
§E2]) By rest-frame day 2797, the Mg Il+Fe II UV62 blend con- 
tains only low- and high- velocity Mg II; Fe II UV62 has vanished. 



typically uncertain by ±0.4 (as compared to ±0.02 for SDSS 
magnitudes). The conversions to SDSS magnitudes were 
made using equation 2 of Monet et al. (2003) except for the 
Palomar Quick V which used Table 7 of Smith et al. (2002). 
In five spectroscopic epochs, the wavelength coverage was 
sufficient to enable synthesis of g — r and r — i colours using 
the SDSS system response curves of Doi et al. (2010) but slit 
losses prevented the synthesis of accurate magnitudes. Only 
in the SDSS spectroscopic epoch could both colours and 
magnitudes be reliably synthesized. Corrections for Galac- 
tic extinction have not been applied to these magnitudes. 

Despite the uncertainties in the magnitudes in Table 2, 
the g — r colours obtained from them provide good evidence 
that the g-band absorption in FBQS J1408±3054 was ap- 
proximately constant in strength for at least 27 rest-frame 
years prior to its weakening sometime after MJD 51612. 
However, we cannot rule out the disappearance and reap- 
pearance of an Fe n BAL structure during that time, espe- 
cially because there is a gap in the historical photometric 
record longer than twice the timescale over which the Fen 
absorption was later observed to vanish. 

In the two SDSS photometric epochs, there was an 18% 
increase in p-band flux relative to the i-band flux between 
MJD 53110 and 53795 (371 rest-frame days). Assuming that 
flux increase was uniform across the g band, we can predict 
that at MJD 53110 (rest-frame day 1759) the Mgn±Fell 
trough had W r = 46 A and the Mgn trough W r = 21 A. 
These points are plotted in Figure|3]without error bars. They 
suggest that the change in W r values during the weakening 
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Table 1. Spectroscopic Observations of FBQS 1408+3054 



Date 


MJD 


Source 


A/AA 


C 




1995-05-23 


49860 


KPNO 4m 


1000 







1996-01-11 


50093 


Lick Kast 


835 


126 


126 


1997-04-10 


50548 


KPNO 2.1m 


1250 


372 


246 


1997-05-16 


50584 


Lick Kast 


835 


392 


20 


1997-06-10 


50609 


Lick Kast 


835 


405 


13 


2000-03-09 


51612 


Keck ESI 


6750 


948 


543 


2005-01-07 


53377 


Keck ESI 


6750 


1903 


955 


2006-03-01 


53795 


SDSS DR7 


1850 


2129 


226 


2009-07-16 


55028 


HET LRS 


1130 


2797 


668 



a Rest-frame days since the first spectroscopic epoch. 

^ Rest-frame days since the previous spectroscopic epoch. 

of the BAL troughs may have been more gradual than rapid, 
and perhaps not monotonic. 

To add to our picture of FBQS J1408+3054, we also 
searched for additional data at other wavelengths (Fig- 
ure SJ. We checked for sensitive X-ray coverage, but none 
is available. FBQS J1408+3054 was detected by 2MASS 
(|Skrutskie et al. 2006|) in 1999 (MJD 51517). With GALEX 
( |Morrissey et al. 20070 , FBQS J 1408+3054 was detected in 
the NUV filter in 2005 (MJD 53523.5) and in both the NUV 
and FUV filters in 2007 (MJD 54224.5). The flux decrement 
between the NUV and FUV filters is consistent with the ex- 
pected strong absorption in the Lyman transitions at the 
quasar redshift. The FUV non-detection on MJD 53523.5 
is consistent with the fainter NUV flux at the same epoch. 
Both GALEX observations occurred after the MJD 53377 
Keck spectrum which shows a strong reduction in Fell and 
Mg II absorption strength, so we cannot determine if the ab- 
sorption at observed NUV wavelengths has weakened like 
that at observed optical wavelengths has. The MJD 53523.5 
NUV flux level is consistent with the MJD 53377 Keck spec- 
trum, and the MJD 54224.5 NUV flux level is consistent with 
a continued decrease in absorption. 



3 BLACK HOLE MASS ESTIMATE 

To estimate the black hole mass of this quasar using the 
FWHM (Full Width at Half Maximum) of its H/3 emission 
line, we must remove the neighboring Fe II emission. 

We subtracted Fe II emission from the MJD 53377 Keck 
ESI spectrum using two templates: Boroson & Green (1992; 
BG92) and Veron-Cetty et al. (2004; VJV04), both con- 
volved with Gaussians of FWHM spaced every 250 km s~ x 
from 1000 to 20,750 km s _1 . First, we subtracted a power- 
law continuum fit made using observed-frame windows 7007- 
7079 A and 10280-10351 A (rest-frame 3791-3830 A and 
5562-5601 A)0Next, each velocity-broadened Fell template 

2 This fit lies above the observed spectrum around 4685 A, sug- 
gesting that FBQS J1408+3054 is one of the rare class of Balmer- 
line BAL quasars (Wang et al. 2008, and references therein). Fell 
subtraction confirms that a significant H/3 trough exists. The 
trough is blueshifted by 11,000 km s _1 , so it does not signifi- 
cantly affect the measurement of the H/3 FWHM. Although the 
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Figure 4. Rest- frame ultraviolet /optical spectrum and rest- 
frame ultraviolet and near-infrared photometry of FBQS 
J1408+3054, plotted as F\ versus rest frame wavelength in A. 
The green points are the 2MASS data from 1999. The blue spec- 
trum is the Keck data from 2005. The red point is the GALEX 
NUV data from 2005. The black points are the GALEX FUV and 
NUV data from 2007. The error bars in the horizontal direction 
for the GALEX points were obtained by taking the difference of 
the effective wavelength and the limits of the passband in both 
directions and averaging. 



was scaled using a progressively increasing normalization 
factor and the \ 2 value between the spectrum and each nor- 
malized template was measured using the rest-frame region 
4440-4647 A. The normalization factor with the lowest x 2 for 
any velocity FWHM was chosen; the best-fit Fe II FWHM 
was 3000 km s _1 in both cases. 

The resulting Fe n-subtracted spectra were used to mea- 
sure the H/3 FWHM. After subtraction of the VJV04 tem- 
plate, we obtained FWHM=4100 km s _1 . After subtraction 
of the BG92 template, we obtained FWHM=3800 km s" 1 . 

To determine the black hole mass from H/3, we use 
the formula Mbh = 5.5FWHM 2 i? Wr /G from Onken et al. 
(2004), where G is the gravitational constant. To determine 
Rtir we use the formula in Bentz et al. (2009): 

log 10 R blr = -2l.3±H + 0.5191^66 log 10 (AL A (5100)) (1) 

where the H/3 broad-line region radius Ru r is in light- 
days, A = 5100 A rest-frame and La (5100) is the 
monochromatic luminosity there. We measure f\ (5100) = 
1.5 x 10 -16 ergs s _1 cm -2 A -1 , which corresponds to 
A£ A (5100) = 2.66 x 10 45 ergs s _1 in a Universe with Ho = 70 
km s _1 Mpc" 1 , Qm = 0.3 and J7a = 0.7. These quantities 
yield R b w = 190 light-days (0.16 parsecs) or 4.9 x 10 12 km. 

We thus find Mbh = 3.4 x 10 9 M e using the VJV04 
template and Mbh = 2.9 x 10 9 M Q using the BG92 tem- 
plate. We adopt Mbh = (3.15 ± 0.35) x 10 9 M , which 
yields a Schwarzschild radius of Rs c h = 9.3 x 10 9 km, so 



H/3 BAL appears to be deeper on MJD 51612 than on MJD 53377, 
we cannot reliably claim a detection of variability. The relative 
flux calibration uncertainties are larger this close to the long- 
wavelength end of the spectrum, and apparent variability in the 
Fe II emission could affect the appearance of the H/3 trough, but 
no H/3 continuum fit or Fell subtraction is possible for the MJD 
51612 spectrum due to its limited wavelength coverage longward 
of H/3. 
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Table 2. Observed and Synthesized Photometry of FBQS 1408+3054 



Date 


MJD 


Source 


Original 


g 


r 




g-r 


r— i 


1 Q^fl OA 1 s 


OOOOiJ 


POSS-I 


0=19.2 


18.6 












OOOOiJ 


POSS-I 


E=17.1 




17.4 




1 2 






45110 


P 1 r\ w~i a v 

JT dlUllldjl 


V=17.5 


18.3 










1 QQ9 C\A OS 


48720 


POSS-II 


j— 1^9 o 


18.8 










1 QQ9 nzL 9R 


48738 


POSS-II 


F=17.3 




17.4 




1 4 




1995-02-23 


49771 


POSS-II 


N=17.0 






17.4 






1996-01-11 


50093 


Lick" 










0.95 


-0.05 


1997-04-10 


50548 


Lick a 










0.78 




1997-05-16 


50584 


Lick a 










1.28 


-0.01 


1997-06-10 


50609 


Lick a 










1.15 


0.05 


1999-12-05 


51517 


2MASS 














2000-03-09 


51612 


Keck° 










0.95 


0.09 


2004-04-13 


53108 


SDSS 


g>r,i 


17.82 


17.40 


17.33 


0.42 


0.07 


2004-04-15 


53110 


SDSS 


g>r,i 


17.84 


17.40 


17.34 


0.44 


0.06 


2005-01-07 


53377 


Keck a 










0.48 


0.07 


2005-06-02 


53524 


GALEX 














2006-03-01 


53795 


SDSS^, 




17.76 


17.47 


17.44 


0.29 


0.03 


2007-05-04 


54225 


GALEX 















Synthesized from spectroscopy; see text. 



that Rur = 530Rsch- Using the mean bolometric correc- 
tion of Richards et al. (2006), we estimate an unremarkable 
Lboi/L e dd = 0.07 ± 0.02 for this object. 



53110 yields gri magnitudes within ±10% of gri magnitudes 
synthesized from SDSS spectroscopy taken on MJD 53795. 



4 DISCUSSION 

The trough variations in FBQS J 1408+3054 are likely due 
to motions of absorbing gas transverse to our line of sight 
(e.g., Hamann et al. 2009). Such motions are seen in numeri- 
cal simulations of line- driven outflows, both two-dimensional 
(Proga, Stone & Kallman 2000; Proga & Kallman 2004) 
and three-dimensional (Kurosawa & Proga 2009). In this 
section we first discuss why a changing ionization explana- 
tion for the trough variations is unlikely. Then we combine 
an estimate of the ultraviolet (UV) continuum source size 
in this object with the timescale over which the absorption 
disappears to estimate the transverse velocity of the absorb- 
ing gas. We then outline a simple BAL outflow model with 
which we constrain the distance of the absorbing gas from 
the black hole and the time since the gas was launched as 
part of the outflow. We close by discussing the implications 
of our results for the connections between BAL subtypes, 
for the average episodic lifetimes of BAL outflows and for 
the link between FR II radio morphology and BAL quasars. 



4.1 Was the Fen Ionized Away? 

One reason it is unlikely that the changes in the BAL troughs 
of FBQS J 1408+3054 are due to a increased ionization of 
the absorbing gas is because the near-UV continuum flux of 
this quasar did not increase significantly while the absorp- 
tion troughs weakened dramatically. The 1997 spectrum pre- 
sented in Becker et al. (2000) and the 2006 SDSS spectrum 
show the same peak flux density near Mgn to within 10%, 
and SDSS photometry of FBQS J1408+3054 taken on MJD 



Of course, ionization is caused not by near-ultraviolet 
photons but by extreme UV and soft X-ray photons with 
A ^ 912 A, and along BAL quasar sightlines such photons 
are thought to be absorbed by a layer of shielding gas (e.g., 
Gallagher et al. 2006). It is possible that a decrease in the 
amount of shielding gas along our line of sight could have 
increased the ionizing flux without increasing the near-UV 
flux that we have observed, resulting in the ionization of the 
Fell absorber to Fein. Increasing ionizing flux would have 
expanded the Fe n ionization front along our line of sight out- 
wards. The timescale for that expansion is the timescale for 
the decrease in the shielding gas column or the recombina- 
tion timescale in the Fe n absorber, whichever is longer. The 
recombination timescale is (cvn e ) _1 , where a=1.02x 10~ 12 
cm" 3 s" 1 is the total recombination coefficient of Fell at 
T=10 4 K (Woods et al. 1981) and n e is the Fell absorber's 
electron density. 



We cannot definitively rule out a change in ionizing 
flux as the explanation for the disappearance of Fell ab- 
sorption in FBQS J1408+3054 because we lack X-ray data 
bracketing the disappearance epoch. We do not believe such 
an explanation is likely, however, because no evidence for 
ionization-dependent absorption variations has been found 
for the general population of BAL quasars (Gibson et al. 
2008) . Nonetheless, if a decreasing shielding gas column was 
responsible for Fe II vanishing from the spectrum of FBQS 
J1408+3054, we can set a lower limit on n e in the Fell ab- 
sorber of between 6 x 10 3 cm -3 and 5 x 10 4 cm -3 correspond- 
ing to the decrease occurring on the maximum or minimum 
observed timescales, respectively. 
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4.2 Constraints on the BAL Structure Size and 
Transverse Velocity 

We refer to the gas that absorbs the quasar's light in Mgn 
and Fell as the BAL structure. Our conclusions about its 
size and location hold regardless of whether the structure 
is a single large cloud, part of a continuous flow tube (e.g., 
Arav et al. 1999), or a collection of smaller subunits (e.g., 
Hall et al. 2007). The column density of the absorption may 
vary over the structure even in the first two cases above 
(de Kool et al. 2002; Arav et al. 2005; Sabra & Hamann 
2005), but for convenience we assume the structure has a 
sharp edge. Note that the entire BAL outflow can be much 
larger in size than the BAL structure we observe along our 
sightline (including extending to greater distances from the 
quasar) and can have a longer lifetime as well. 

Given the size of the continuum source, the timescale 
over which absorption changes from strong to weak (or vice 
versa) constrains the velocity and location of the BAL struc- 
ture. Furthermore, the length of time the quasar was seen 
as a BAL quasar constrains the size of the BAL structure 
in one transverse direction. For our discussion, we take the 
change in the spectrum of FBQS J1408+3054 that began 
between 2000 and 2005 to be a change from strong to weak 
absorption. 

The minimum timescale over which the Fell troughs 
disappeared is the time between the 2005 and 2006 spec- 
tra: At m i„=226 rest-frame days. The maximum allowed 
timescale is the time between the 2000 and 2009 spectra: 
At max =1849 rest-frame days. We adopt an intermediate 
timescale of At — 946 rest-frame days as the time required 
for the Mg Il+Fe II rest-frame equivalent width to decrease 
from its level in 2000 to its level in 2009 at the rate of change 
observed between 2005 and 2006. 

This value of At = 946 rest-frame days is an estimate of 
the crossing time of the quasar emission region by the trail- 
ing edge of the BAL structure^ As the quasar was seen as 
an FeLoBAL quasar of roughly constant W r for at least 948 
rest-frame days prior to the 2000 epoch, the time required 
for the entire BAL structure observed in FBQS J1408+3054 
since its discovery to cross the quasar emission region was 
at least t croas = 2At. In other words, the entire BAL struc- 
ture took at least twice as long to cross the quasar emission 
line region as the edge of the BAL structure didQ Twice 
the 2700 A continuum emission region size D2700 is thus an 
approximate lower limit on the transverse size of the BAL 
structure in the direction of its projected motion on the sky. 
The maximum depth of the Fe 11 absorption reveals that the 

3 We have assumed a sharp-edged structure, but the results apply 
generally. Because the Fe II absorption goes from its maximum to 
zero in ~946 rest-frame days, a diffuse structure in which TVjf 
gradually changes at the edges would have to move across the 
face of FBQS J1408+3054 faster than a sharp-edged structure. 
Therefore, the minimum transverse velocity derived for a sharp- 
edged structure is a lower limit to the velocity of any structure. 

4 The lack of strong variability in the historical photometry of 
FBQS Jf408+3054 discussed in § \T2\ suggests that the BAL 
structure actually had a crossing time of ;J 27 rest-frame years, 
~10 times longer than the crossing time of its edge. However, 
as we lack spectroscopy or regularly sampled photometry in pre- 
discovery epochs, we quote the shorter timescale to be conserva- 
tive. 



BAL structure covered at least ~75% of the region that pro- 
duces the continuum emission at 2700 A (and possibly part 
of the Fe 11-emitting region, but not all; see below) . Because 
this covering fraction was roughly constant for several rest- 
frame years, we take it to be the extent of the BAL structure 
in the direction on the sky orthogonal to the structure's mo- 
tion on the sky. 

We use a Shakura & Sunyaev (1973) thin accretion disk 
model to estimate -D2700, the diameter of the disk within 
which 90% of the 2700 A continuum is emitted, for a Mbh = 
3.15 x 10 9 M Q black hole. We find D2700 = 4.28 x 10 11 km 
(AQRsch)- The projected dimensions of the BAL structure 
are therefore Si 3.2 x 10 11 km or 0.01 pc wisde by Si 8.6 x 
10 11 km or 0.028 pc (and likely > 8.6 x 10 12 km or 0.28 pc) 
long in its projected direction of motion. 

Using At and D2700, we can constrain the velocity of 
the BAL structure perpendicular to our line of sight: v± = 
D2700 / At = 5200l26oo k m s_1 > where the range comes from 
the constraints on At. The uncertainty on the size of the 
region crossed by the flow can change these limits somewhat 
(see the end of the next section). Nonetheless, the transverse 
velocity is constrained to be approximately 2600 km s" 1 
^ ki ^ 22,000 km s _1 ); for comparison, the line-of-sight 
velocity of peak absorption is a blueshift of 12,000 km s _1 . 

We can rule out the possibility that the BAL structure 
completely crossed the larger Fe 11-emitting region during 
the time span At. If we use 27?Feii instead of D2700 as the 
size of the region crossed by the BAL structure, v±_ will in- 
crease. As the Fell FWHM in this object is smaller than the 
H/3 FWHM, Fe 11 emission is likely to come from a larger av- 
erage radius than H/3. Assuming R oc FWHM -2 , we obtain 
R Fe n = l.lARbir = 8.5 x 10 12 km = 910i? Sch . In that case, 
we find v t , F eii = 2R FcU /At = 206, 000 km s" 1 . Such a high 
velocity implies a small distance from the quasar, which is 
not self-consistent. A BAL structure located at any distance 
d-BAL < -Rpdi would not be able to cover the full extent of 
the Fell emission region. 

4.3 A Model for Constraining the BAL Structure 
Location 

When the radial velocity of absorbing gas around a quasar 
is not known, its distance from the black hole is usually es- 
timated by assuming it is in a circular orbit (e.g., Risaliti et 
al. 2007). However, in our case the radial velocity of the BAL 
gas is 12,000 km s _1 and we cannot self-consistently assume 
a circular orbit. To estimate the distance of the BAL struc- 
ture from the black hole, cLbal, some assumptions about the 
velocity field of the BAL wind must be made (see, e.g., § 6.2 
of Misawa et al. 2005). 

Given the large blueshift of the BAL structure in FBQS 
J1408+3054, the outflow was likely launched at a radius 
where the escape velocity is comparable to the observed 
blueshift (e.g., Murray et al. 1995; Proga, Stone & Drew 
1998; Proga 1999; Everett 2005; Zhang & Thompson 2010). 
A natural source for such gas is a wind from the accretion 
disk around the quasar, and that is the scenario we adopt. 

We assume that our line of sight to the quasar lies at a 
latitude A above the accretion disk and that the accretion 
disk has azimuthal symmetry (Figure [SJ). The BAL wind 
crosses our line of sight to the continuum emission region 
at a distance cLbal from the black hole, at cylindrical co- 
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ordinates tbal = (IbalcosA. and zbal = Abal sin A. We 
parametrize the velocity field of the BAL wind at this point 
as an azimuthal velocity V<f, and a poloidal velocity V w ori- 
ented at an angle •& above the accretion disk[f| The BAL 
velocity components parallel and perpendicular to our line 
of sight are 



Wy = Vw cos(i? — A) 

v_l = ^V| + V2sin 2 (0-A). 



(2) 
(3) 



We know Vn from spectroscopy, and the BAL variability con- 
strains V ± , min < Vj_ < Dl, m „. 

To link these velocities to cLbal, we first assume that 
the BAL outflow started out in a circular orbit at radius 
ro in the accretion disk and has not swept up a significant 
amount of mass since then. We assume conservation of an- 
gular momentum in the <j> direction (neglecting momentum 
carried by magnetic fields) so: 



V<t> = V GMbh ro I tbal 



(4) 



where tbal is the current distance of the BAL structure 
from the axis of symmetry of the accretion disk. 

Next, we use the results of Murray et al. (1995, fig. 
3) and Murray & Chiang (1998, equation 2) that the final 
poloidal velocity of a radiatively driv en outflow launched 
from radius r is Voa = (3.6 ± 1.1)^J GMbh /ro and that 
Voo is reached within a radial distance of a few times ro. 
Therefore, we assume 



V W =V X = (3.6 ± 1.1) y/ GMbh /ro 
and write 



«H = (3.6 ± 1.1) J ° MbH cos(i9 - A) 



ro 



(5) 



(6) 



v ± = 



GMbhto 



(3.6 ± 1.1)3 GM BH s . n2 ^ _ ^ ^ 



ro 



These equations express Vu and v±_ in terms of physical 
parameters in our outflow model, in principle enabling us 
to constrain those parameters and place limits on cLbal = 
rBAh/ cos A. 

In the case of known A and ro can be obtained from 
equation [6] A lower limit on v± then yields an upper limit 
on rsAL via equation!?] an d an upper limit on cLbal follows. 
Similarly, a lower limit on cLbal follows from an upper limit 
on v±. 

Note that if A > #, the streamlines must have bent to 
reach the location (tbal, zbal) with velocity V w at angle 
The streamlines therefore must also have crossed our line 
of sight at another position closer to the black hole (Fig- 
ure |SJ right-hand panel), possibly producing another BAL 
trough in the spectrum. Another trough is guaranteed to be 
produced only if the BAL outflow is continuous back to the 
launch radius and extends over all azimuthal angles. 

We make one final refinement to the above analysis. We 
have assumed that continuum emission region is intrinsically 
circular with diameter D2700, projected on the sky to an 



5 The magnitude of the three-dimensional velocity v of the BAL 
wind at this point is given by v = ^/V^ + V% = \J v \ + ■ 



ellipse of major axis -D2700- The transverse velocity limits 
were computed assuming a distance D2700 was crossed by 
the BAL flow, but in fact the trailing edge of the BAL flow 
almost always crosses a distance D± < D2700 (Figure [6]). 
D± is defined as the distance spanned by the projection of 
the ellipse onto an axis parallel to v± ; only when A = # does 
D± = D2700- For a given $ and A, once we calculate ro from 
equation [6] we calculate the direction of v± on the sky and 
the distance D± along that direction, then recalculate 



V±, r , 



(8) 



for comparison to equation [7J and iterate to convergence. 



4.4 Constraints on the BAL Structure Location 

Unfortunately, we only have statistical constraints on A and 
We must therefore consider what constraints we can put 
on ro and r bal given reasonable distributions of our viewing 
angle A and the wind velocity's poloidal angle above the disk 
plane 

We assume 0° < A < 73° from the size of this FR II 
quasar's radio lobes relative to the maximum size known for 
such lobes (Gregg et al. 2006). (Even without identifying this 
quasar as an FR II quasar, we can assume A ^ 80° because 
Hovatta et al. (2009) find that objects viewed at A > 80° are 
seen as blazars.) We assume 5° ^ ■& ^ 85° for the outflow 
angle because we have no real constraints on this parameter. 
We sample both angles using uniform distributions in solid 
angle. 

We compute lower and upper limits on cLbal for each 
combination of A and 1? spanning the above limits, and 
take the averages. Combinations for which equation [7] yields 
the square root of a negative number as an upper limit 
on v± are ruled out in this model and do not place any 
constraints on cLbal- Also ruled out are combinations that 
yield r B AL,max < r for $ < A or r B AL,max < r /(l - 
tan A/ tan •&) for ■& > A. 

There are two sources of uncertainty in the lower and 
upper limits we find. The first is the scatter in results be- 
tween different #,A combinations. This uncertainty can be 
reduced with better knowledge of quasar viewing angle dis- 
tribution^] or FeLoBAL wind outflow angles. However, even 
without such knowledge, the ratio of upper to lower limits on 
the BAL distance is well constrained; the values of such lim- 
its vary considerably with 1? and A, but their ratio does not. 
The second is the uncertainty in inferring ro from the factor 
of 3.6±1.1 in equationJS] The larger this factor, the farther 
out the wind can be launched and still end up with the same 
observed velocity, for any i9,A combination. Thus, the BAL 
distance increases with this factor. This uncertainty will re- 
main unless studies of disk winds can refine the distribution 
of their asymptotic wind velocities. 

Including the above uncertainties, we find average lower 



and upper limits on (Ibal of 5800 



+9400 



< di 



■JR. 



6 For example, we have assumed a flat probability distribution of 
the viewing angle A in FBQS J1408+3054, 0°<A<73°. A more 
realistic probability distribution could be constructed by compar- 
ing the observed separation between its radio lobes to an inferred 
distribution of true lobe separations, or by using other methods 
of inferring radio source orientations (Richards et al. 2001). 
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Figure 5. Side view of a BAL structure (open circle) obscuring the our line of sight (L.O.S.) to a black hole and accretion disk (filled 
circle at the origin and thick horizontal line segment). Our line of sight is at an angle A above the disk. The BAL structure at location 
{tbal, zbal) moves at velocity V w at angle i3 above the disk in the r — z plane and at azimuthal velocity normal to the r — z plane. 
The dashed line segments show a plausible upper limit on the projected distance the BAL wind has traveled from its launch radius. Two 
geometries that yield equal and opposite i? — A arc shown: on the left, = 50° and A = 35°; on the right, i? = 35° and A = 50°. 
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Figure 6. View from along our line of sight of a BAL structure (large rectangle extending off the plot) partially obscuring the continuum 
emission region around the black hole. The continuum emission region, indicated here with horizontal lines, is intrinsically a circular disk 
but is projected to an ellipse. The BAL structure moves with velocity vj_ , which can be broken down into two components: the azimuthal 
velocity and the component of the wind velocity in the plane of the sky V w sin(i9 — A). The same two geometries as in the previous 
figure arc shown: on the left, i? = 50° and A = 35°; on the right, i? = 35° and A = 50°, in which the emission region appears more 
circular. The width of the continuum emission region projected along the v± direction, D±, is indicated by the dashed line segment; D± 
is slightly larger in the geometry shown on the right. Although the values of 1? and A are unknown in FBQS J1408+3054, these figures 
are otherwise representative of FBQS J1408+3054, where we observed the disappearance from along our line of sight of an Fell BAL 
structure with a minimimum size larger along v± than perpendicular to it and which even at maximum absorption likely covered only 
~ 75% of the continuum source. 



Variable BAL Quasar FBQS J 1408+3054 11 



Table 3. BAL Distance Limits, Averaged Over i? and A 



tanw < 



Tbal tan A 



tbal 



Hi 



(l-smfi)+r B AL tanA.(lO) 



\/GM BH /r 



d-BAL, 



-) (Ratio) 



2.5 
.3.6 

4.7 



2goo +l800 

5800^°° 460001^000 
990011%° 7100011°*°°° 



290001*?°°° 9.1±0.2 
' ^ 7.0±0.2 



i.3±0.2 



46000l^gg° (90% confidence ranges). While these limits 
are each uncertain by more than an order of magnitude, the 
ratio of upper to lower limits has a much more well defined 
value of ~7.5±1.5. 

The large uncertainties are in part inherent to our lim- 
ited knowledge of the orientation of the outflow, but are 
also influenced by the fact that the variability timescale 
(and thus the transverse velocity) in this object is only con- 
strained to within a factor of eight. Better limits on the 
variability timescale would yield tighter constraints on the 
BAL location and restrict the A combinations which could 
match the observations. 

For reference, in Table 3 we list the average lower and 
upper limits on the distance of the BAL structure in FBQS 
J1408+3054 as a function of the assumed ratio between the 
BAL's final and initial poloidal velocities. The 90% confi- 
dence limits in the Table come only from the scatter in the 
limits for different §, A combinations. 

These limits on ds al place the BAL structure in FBQS 
J1408+3054 between 11 and 88 times farther from the 
quasar than the H/3 broad-line region. Physically, these lim- 
its correspond to a distance between 1.7 and 14 pc from the 
black hole. 

Combinations of $ and A that match the observations 
have an average value of $ — A ~ —16°, with a range of 
-60° < ■& - A < 41° out of the allowed range of -68° < 
i? — A < 85°. The angle •d — A and the magnitudes of V w 
and V$, constrain the angle if) between the line of sight and 
the velocity vector of the wind. For combinations of $ and A 
that match the observations we find an average if) ~ 24° ±2°, 
with a range of 8.9° < tp < 60° out of the allowed range of 
1.2° < V < 61.4°. 



4.5 Constraints on the BAL Structure Age 

The age of an observed BAL structure (how long it has ex- 
isted) is a lower limit on its lifetime (how long it will exist) 
and on the lifetime of the entire BAL outflow. A BAL struc- 
ture's age can be constrained by considering how long gas in 
the structure took to reach its current location after being 
launched. This constraint applies regardless of whether the 
structure is the result of a single ejection event or part of a 
continuous flow. 

To constrain the age of the BAL structure, we trace the 
wind at angle •d from (tbal, zbal) back to r = ro or z — 0, 
whichever comes first, and then back to (ro,0); see Figure 
[5] Such a path yields the longest plausible distance s for the 
BAL structure to have covered to reach its current position. 
The formulae for s are: 

, r S AL tanA tanA 

tanw ^ : s = tbal— — r(l— cost))+rBAL— ro(9) 

tbal — ro sin § 



tbal — ro COS l3 

The age of the BAL structure is then simply t age — s/V m . 
This value will be a slight underestimate because the wind 
was launched with v <V W , but the underestimation will be 
small for s > ro. 

For combinations of i? and A that match the obser- 
vations of FBQS J1408+3054, we find average BAL struc- 
ture ages of tage = 901^ years for cLbal = 5800 R ac h and 



t aqe = 1100: 



years for cIbal = 46000 R ac h- These lower 



limits on the BAL outflow lifetime in FBQS J 1408+3054 are 
consistent with other, independent determinations (see § 12.21 
and §H7b. 



4.6 Implications for Links Between BAL Subtypes 

Aside from their defining broad absorption at UV wave- 
lengths and accompanying X-ray absorption (e.g., Gallagher 
et al. 2006), HiBAL quasars as a population do not appear 
distinct from non-BAL quasars in terms of their multiwave- 
length properties (Lewis et al. 2003; Willott et al. 2003; Prid- 
dey et al. 2007; Gallagher et al. 2007; Shen et al. 2008). The 
only exception is the greater reddening of HiBAL quasars 
relative to non-BAL quasars (Brotherton et al. 2001; Re- 
ichard et al. 2003; Trump et al. 2006; Gibson et al. 2009). 

Numerous studies have suggested that LoBAL and 
FeLoBAL quasars may have some optical and mid-infrared 
(mid-IR) spectral properties distinct from HiBAL and non- 
BAL quasars, such as weaker [O ill] emission, stronger Fell 
emission and continuum polarization in the optical and 
ultraviolet, and higher mid-IR luminosities and different 
mid-IR spectra (for LoBALs and the low-redshift FeLoBAL 
quasar Mrk 231, see Low et al. 1989, Weymann et al. 1991, 
Boroson & Meyers 1992 and Canalizo & Stockton 2001; for 
both LoBALs and FeLoBALs, see Schmidt & Hines 1999 and 
Urrutia et al. 2009; for FeLoBALs, see Farrah et al. 2007, 
2010). On the other hand, some studies have found little 
evidence that LoBAL or FeLoBAL quasars have properties 
significantly different from those of HiBAL quasars (for po- 
larization properties, Ogle et al. 1999; for mid-IR properties, 
Lewis et al. 2003; Willott et al. 2003; Gallagher et al. 2007). 

The transformation of FBQS J 1408+3054 from an 
FeLoBAL to a LoBAL quasar along our line of sight casts 
these studies in a somewhat different light. Suppose that at 
least some FeLoBAL quasars have properties distinct from 
HiBAL and LoBAL quasars in some wavelength range where 
the timescale for variability is longer than in the optical (e.g., 
in the mid-IR; Farrah et al. 2007, 2010). Then the transfor- 
mation of FBQS J 1408+3054 from an FeLoBAL to a LoBAL 
requires that some, but not necessarily all, LoBAL quasars 
must have the same distinct properties. In this scenario, the 
reason why some but not all LoBAL quasars have unusually 
high mid-IR luminosities as compared to HiBAL quasars has 
two underlying causes. First, at least some FeLoBAL quasars 
have such high mid-IR luminosities. Second, some objects 
that are seen as LoBAL quasars along our line of sight could 
be classified as FeLoBAL quasars when seen along some lines 
of sight, while some LoBAL quasars would not be seen as 
FeLoBAL quasars along any line of sight. 

If FeLoBAL quasars as a population do have unusual 
multiwavelength properties, they are unlikely to be "nor- 
mal" BAL quasars seen at specific viewing angles. They 
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could instead be quasars emerging from a shroud of gas 
and dust (Voit et al. 1993), for example as transition ob- 
jects between ultraluminous infrared galaxies (ULIRGs) and 
quasars (Farrah et al. 2010), or quasars with a specific 
disk wind geometry (Richards 2006). Our study of FBQS 
J1408+3054 does not provide support for the transition ob- 
ject scenario because the FeLoBAL gas in this object moves 
at very high velocities and is found only 1.7 to 14 pc from 
the black hole. If FeLoBALs are to be transition objects 
between ULIRGs and quasars, they must drive out gas on 
kiloparsec and not just parsec scales. That said, our study 
does not rule out this scenario either. Distance constraints 
on more FeLoBAL outflows would be one way to deter- 
mine how common transition objects are among FeLoBAL 
quasars. Another possible test of the transition object sce- 
nario is to determine if FeLoBAL quasar host galaxies differ 
from the host galaxies of other quasars, due to more recent 
interactions or mergers. 

However, such differences are not necessarily required in 
the original 'clearing shroud' scenario of Voit et al. (1993). 
That work postulated that LoBAL quasars could be quasars 
emerging from parsec-scale shrouds of dusty gas, and our ob- 
servations are entirely consistent with such a scenario. The 
existence of dusty gas on such small scales need not depend 
on the kiloparsec-scale structure of the quasar host galaxy, 
so FeLoBAL quasar host galaxies could be unremarkable in 
this model. 

In the alternate scenario where FeLoBALs are quasars 
with specific disk wind geometries, they must be differ- 
ent from other quasars in at least one parameter driving 
accretion disk structure. For example, while the distribu- 
tion of Eddington ratios is known not to differ between 
non-BAL and HiBAL quasars (Ganguly et al. 2007), the 
FeLoBAL quasar Eddington ratio distribution has not been 
studied; however, FBQS J1408+3054 has an unremarkable 
L bol /L Edd = 0.07 (§|3. 

It is known that greater average reddening is found 
among ever rarer BAL subtypes. That is, as mentioned 
above HiBAL quasars are on average more reddened than 
non-BAL quasars, LoBAL quasars are more reddened than 
HiBAL quasars (Sprayberry & Foltz 1992), and FeLoBAL 
quasars are more reddened than LoBAL quasars (Re- 
ichard et al. 2003). If variability of the sort seen in FBQS 
J1408+3054 is common, then we predict that quasars of one 
subtype which are more reddened than average for that sub- 
type will be more likely to develop BAL troughs of the next 
subtype. (Note that when its SDSS spectrum was taken, 
FBQS J1408+3054 itself had A(g - i) = 0.33, bluer than 
the average A(g — i) ~ 0.5 for the LoBAL quasars from Re- 
ichard et al. (2003). A(g — i) is the g — i colour of a quasar 
minus the average g — i for quasars at the same redshift.) 

4.7 Average Episodic Lifetime of BAL Outflows 

We can use multiple-epoch spectroscopy of BAL (and non- 
BAL) quasars to constrain the average episodic lifetime 
of BAL outflows by expanding upon §4.2 of Gibson et al. 
(2008). We take all time intervals in the quasar rest frame. 

We define an episodic lifetime as the time over which a 
quasar's spectrum exhibits a BAL outflow, using a predeter- 
mined definition of what constitutes a BAL trough. (Note 
that by this definition, BAL episodes do not overlap in time.) 



The appearance and disappearance of a BAL in a spectrum 
thus define the start and end of a BAL episode, regardless of 
whether they reflect the creation and destruction of a BAL 
outflow or merely the passage of one or more structures in 
a BAL outflow across our line of sight. The episodic BAL 
lifetime is what we can observe, but it is worth remembering 
that the lifetime of the entire BAL outflow must be equal 
to or longer than the episodic lifetime while the lifetime of 
a structure in a BAL outflow can be longer or shorter than 
the episodic lifetime (§[43}. 

Suppose that quasars are sufficiently luminous to be 
included in a sample for an average time tQ (differential se- 
lection effects would admittedly have to be considered unless 
bolometric luminosity selection could be used) that a frac- 
tion fhos of those quasars exhibit an average of Nb BAL 
episodes along our line of sight during their lifetimes (mean- 
ing that a fraction 1 — jhos never do so), and that those 
BAL episodes have an average episodic lifetime ts- If we 
consider ts and tQ to be fixed, although they are actually 
the means of distributions, then the average instantaneous 
fraction of quasars showing BALs along our line of sight is 



fBAL = fLOsNstB/tQ. 



(11) 



Gibson et al. (2008) assumed }los = 1; we relax that as- 
sumption to /ios ~ 0.2, where the lower limit is set by the 
observed BAL fraction fobs — 0.2 (e.g., Knigge et al. 2008). 
Note that equationQTJdoes n °t depend on how BAL episodes 
are spaced over the quasar's lifetime, as long as the episodic 
BAL lifetime ts does not change systematically over the 
quasar's lifetime. 

Now consider two observations of a BAL quasar sepa- 
rated by some (rest-frame) time interval. 

The probability p~ of observing a BAL disappear from 
a BAL quasar over a time interval t~ bs < ts is 

P~ = tobs/tB (12) 

if we also have t~ bs < t gap , where t gap = (tQ/N B ) — t B is the 
average time between BAL episodes. 

The probability p + of observing a BAL appear in a non- 
BAL quasar over a time interval t^ bs < is is conceptually 
similar. It is the fraction of all non-BAL quasars which ever 
appear as a BAL quasar ( ^ L °^fJf^ L ) times the number 
of BAL episodes that start during such quasars' observable 
lifetimes (~ Nb) times the probability of a BAL episode 
starting during time t^ bg (~ t^ bs /tQ). To write a more ex- 
act expression for p + , we must account for the time non- 
BAL quasars which can appear as BAL quasars spend as 
BAL quasars, and for the fraction fo of such quasars which 
have BALs at the start of their observable lifetimes. The full 
derivation is given in the Appendix; the resulting expression 
for p + is: 



fBAL t a 



A 

N B 



(13) 



1 — fBAL ts 

Neglecting the correction term in brackets (unity for fo~0 
or large Nb), if all quasars appear as BAL quasars at some 



point (/los = 1) then we have /sal = 



(equation 



lll[) and we recover p 



That expression is the 



l-N B t B /t Q ■ 

Gibson et al. (2008) approximation p + = N E t^ bs /tQ with 
the quasar lifetime tQ replaced by the time a quasar spends 
as a non-BAL quasar. 
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In practice, by constraining p~ we constrain tB, 
and that constraint plus knowledge of /bal constrains 
fLOsNB/tQ via equation 1111 and }o/Nb via equation 1131 
Limits on tQ can be invoked to constrain JlosNb, but only 
in certain cases might it be possible to constrain fo, Nb 
or /los separately by monitoring the appearance or dis- 
appearance of BAL outflows (see below). The assumptions 
underlying this analysis are that fixed tg and ts are rea- 
sonable approximations and that no systematic change in 
tB occurs over the quasar's lifetime (made to derive Jbal), 
that the timescale t~ bs over which a BAL outflow is observed 
to disappear obeys t~ bs < ts and t ~ bs < t gap (made to derive 
p~) and that the timescale t~ bs over which a BAL outflow 
is observed to appear in a non-BAL quasar obeys t^ b < tB 
(made to derive p + ). Future analyses of BAL episodic life- 
times could relax these assumptions as well as consider BAL 
episodes that overlap in time, whether at the same velocities 
or at different ones. 

Meanwhile, it is instructive to calculate the constraints 
obtained by considering the BALs in FBQS J1408+3054 
to have disappeared, even though the Mgn BAL has not 
quite done so. Those constraints are presented in the follow- 
ing paragraphs, along with more conservative constraints (in 
parentheses) obtained by considering FBQS J 1408+3054 as 
a case of extreme BAL variability, but not of BAL disap- 
pearance. 

Of 13 LoBALs and FeLoBALs discovered in the FIRST 
Bright Quasar Survey (Becker et al. 2000, 2001) and with 
later spectra from the SDSS, only FBQS J1408+3054 shows 
dramatic changes in the absorption covering fraction and 
equivalent width. FBQS J 1408+3054 is still the only such 
strongly variable BAL quasar when the parent sample is 
expanded to include the 24 FBQS HiBAL quasars with later 
spectra from the SDSS. Using the binomial distribution, this 
result translates to a constraint of p~ = 0.0271° 024 (p~ < 
0.060); all limits and ranges given in this section are 90% 
confidence. We have also inspected 156 FBQS non-LoBALs 
with later spectral coverage of the Mgn region from the 
SDSS. None of those quasars exhibited Mgn BALs in their 
SDSS spectra, yielding a constraint of p + < 0.0147. 

Using that combined p~, we constrain individual BAL 
lifetimes to be 45 < ts < 980 years (ts > 60 years), as com- 
pared to the Gibson et al. (2008) result of ts > 43 years. 
This range of lifetimes is broadly consistent with indepen- 
dently derived age limits on the BAL structure in FBQS 
J1408+3054 (§[22]and § [43J), Taking f B AL = 0.2, the com- 
bined constraint on p + yields the constraint fo/Ns > 0.5 
for tB = 45 years but no constraint (fo/Ns > 0) for any 
tB > 97 years. No constraint on fo/Ns is possible if we 
consider FBQS J1408+3054 to be a case of extreme BAL 
variability rather than of BAL disappearance. 

Note that placing a sufficiently strong upper limit on 
the episodic BAL lifetime ts has the potential to place inter- 
esting constraints on fo/Ns- For example, suppose we have 
a constraint ts < 86 years from a p~ = 0.0471q'q^2 mea- 
surement of 14 BAL troughs in a sample of 300 disappear- 
ing over 3 years (a rate within the current 90% confidence 
limits). Then the current constraint of p + < 0.0085 would 
require fo/Ns > §. That would mean Nb ^ 9, as fo ^ 1 
by definition. In turn, that would constrain tg < 3800f los 
years, or tQ ^ 3800 years as /los ^ 1 by definition as well. 
(Such a lifetime is a factor of three smaller than the min- 



imum episodic quasar lifetime suggested by the proximity 
effect; see Martini 2004.) In general, a limit on tQ can be set 
when p~ is significantly larger than p + , as follows. If BALs 
are seen to disappear more frequently than they appear, 
then maintaining a constant fBAL requires that a signifi- 
cant fraction of quasars enter the sample as BAL quasars 
and that Nb not be too large, so that eliminating one BAL 
appearance per quasar has an appreciable impact. Once up- 
per limits are available for Nb and tB, an upper limit to 
tQ follows because the quasar lifetime cannot be arbitrarily 
larger than NstB while still maintaining the observed fBAL- 



4.8 Implications for the FR II-BAL Link 

FBQS J 1408+3054 was identified as likely having an FR II 
radio morphology by Gregg et al. (2006), who found that 
BAL quasars with FR II morphologies are roughly a fac- 
tor of ten less common than would be expected if FR II 
morphologies and BAL absorption were independent phe- 
nomena. FBQS J1408+3054 qualifies as a radio-loud FR II 
quasar if the radio sources identified as its FR II lobes are 
in fact associated with it, or as a radio-intermediate, core- 
dominated quasar if not (Gregg et al. 2006). 

The weakening of the Mgn absorption in FBQS 
J1408+3054 slightly reduces the significance of the anticor- 
relation between BAL strength (as measured by the balnic- 
ity index) and radio-loudness parameter R* for FR II BAL 
quasars found by Gregg, et al. (2006). However, it remains 
true that FR II BAL quasars define the upper envelope of 
that anticorrelation among all radio-selected BAL quasars. 

The implications of the FR II-BAL link for BAL quasar 
variability depend on the reason(s) for the FR II-BAL link. 

There could be some (magneto) hydrodynamic mecha- 
nism which disfavors the contemporaneous production of 
strong BAL winds and strong radio jets that produce FR II 
morphologies. Such a mechanism could explain why HiBAL 
and LoBAL quasars are underrepresented among quasars 
with logi?* > 2 (White et al. 2007). Any such mechanism 
might also have to explain why LoBAL quasars are over- 
represented among quasars with 1 < log R* < 2 (White et 
al. 2007); however, no correction was made in that analysis 
for the substantially greater reddening of LoBAL quasars 
as compared to HiBAL quasars (Reichard et al. 2003). In 
any case, the transience and variability of BAL outflows as 
a function of R* in this model would depend on the column 
density variations of the outflows as a function of R* . 

Alternatively, what if the anticoincidence of BAL winds 
and FR II radio morphologies occurs because the former 
tend to occur early in a quasar's lifetime and the latter can 
occur only when a BAL shroud has been sufficiently reduced 
in column density by the quasar's radiation pressure? In such 
an 'evolutionary' model, quasars with larger lobe separations 
should on average be older and nearer the ends of their BAL 
phases, resulting in greater transience and variability of BAL 
outflows along a given sightline in such objects. 

Observations of BAL variability and transience as a 
function of lobe separation and R* are needed to test the 
evolutionary model and provide constraints any dynamical 
model of BAL outflows must match. 
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5 CONCLUSION 

The quasar FBQS J1408+3054 was discovered as a spectac- 
ular FeLoBAL quasar, and historical data indicate it had 
been one for ^ 20 rest-frame years prior to its discovery, 
but it is now only a modestly absorbed LoBAL quasar. Its 
Fell absorption outflowing at 12,000 km s _1 disappeared 
over a rest-frame time span of between 0.6 and 5 years, and 
its Mg II equivalent width decreased by a factor of two over 
the same time period. This variability was likely caused by 
a BAL structure moving out of our line of sight to the ultra- 
violet continuum emitting region of the quasar's accretion 
disk, indicating a transverse velocity between 2600 km s _1 
and 22,000 km s~\ 

In the context of a disk wind model, we connect the 
observed radial velocity and variability timescale of the Fe II 
BAL structure to its velocity vector and distance from the 
quasar. We have studied what constraints can be placed on 
those quantities given the unknown inclination of the accre- 
tion disk and the unknown angle between the wind veloc- 
ity and the disk plane. We estimate that the BAL structure 
which moved out of our line of sight in FBQS J1408+3054 is 
located between 5800 Rsch and 46,000 Rsch from the black 
hole and that the velocity vector of the wind is oriented at 
~24° to our sightlinc. 

There could be other BAL structures at different dis- 
tances from FBQS J1408+3054, but the close correspon- 
dence of the outflow velocity of Mgll with that of Fell 
strongly suggests that the Mgll absorbing structure is lo- 
cated at the same distance as that of Fell. Continued mon- 
itoring of FBQS J 1408+3054 will be useful to constrain the 
size of the Mg II structure, should it too move out of our line 
of sight, and to constrain the time required for another Fell 
absorption structure to move into our line of sight. 

Lastly, we have worked out in greater detail than hereto- 
fore how multiple-epoch spectroscopy of BAL and non-BAL 
quasars can be used to constrain the average lifetime of BAL 
episodes, ts, and potentially the average number of BAL 
episodes per quasar. At the moment, the 90% confidence 
limit on ts is >60 rest-frame years. Future observations of 
BAL quasar trough variability should greatly improve upon 
that constraint. 
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APPENDIX A: THE PROBABILITY OF 
OBSERVING A BAL OUTFLOW APPEAR 

Equation [TJ] is derived as follows. Consider a quasar of ob- 
servable lifetime tq which has Nb BAL episodes during that 
lifetime, each of duration ts (thus the BAL episodes do not 
overlap). We can observe the quasar as a non-BAL for a total 
time tq — NstB - Each BAL phase is preceded by a non-BAL 
phase (except in quasars which start their observable life- 
times as BALs, for which we make a correction later). Thus, 
^ B ^tbs 16 th e total amount of time which is t^ b later than 
any non-BAL phase and in which the quasar is a BAL. The 
probability of observing a non-BAL as a BAL some time 
^tbs later is the ratio of that amount of time to the total 
amount of time which is t^ bs later than any non-BAL phase: 
tq — NstB — xttbs — ~ NstB- (The x factor is zero for 
quasars which end their observable lifetimes as BALs, and 
unity otherwise. Its average value x depends on the relative 
frequency of those occurrences, about which we do not wish 
to make any assumptions (though we note that for randomly 
occurring BAL episodes, x — 1 — NstB /tg)- Fortunately, our 
assumption that t^ bg <ts makes it reasonable to neglect the 
term involving x.) Thus we write p + as: 



+ _ /los - Sbal N B t+ bs /t Q 



1 - Sbal 1 - N B t B /t Q 



(Al) 



Replacing f LOS with %^ 
factor r 2 -: 



(equation [TT]) and inserting a 



t obs fBAL N B t B 
1 — fBAL tQ tB 



N B t B 



1 - N B t B /tq 



Multiplying 



Nb*i 



p = 



ttbsfBAL/ti 



fE 



through the first term in brackets: 



[1 - N B t B /tQ} 



1 - N B t B /t Q 



(A2) 



(A3) 



The terms in brackets cancel, leaving only the first term. 

Finally, consider a quasar which starts its observable 
lifetime as a BAL quasar. It will have only Nb — 1 BAL 
episodes which could be observed to start. If a fraction fo of 
quasars which exhibit BALs at some point start out as BAL 
quasars, then we have 



P = 



KbsfBAL/tB f (N B - 1) , t+ b J B AL/t B 



fBAL 



N B 



- T" , - (l-/o)ffi(A4) 
1 — JBAL Nb 



tt bs fBAL/tB f N B - fo +N B - f N B 



p 1 f 

t — JBAL 

yielding equation 1131 at last. 



(A5) 
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